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ABSTRACT 

We model the emission lines generated in the photoionised debris of a tidally disrupted 
horizontal branch star. We find that at late times, the brightest optical emission lines 
are [N n] AA6548, 6583 and [O in] AA4959, 5007. Models of a red clump horizontal 
branch star undergoing mild disruption by a massive (50 - 100 Mq) black hole yield 
an emission line spectrum that is in good agreement with that observed in the NGC 
1399 globular cluster hosting the ultraluminous X-ray source CXOJ033831.8 - 352604. 
We make predictions for the UV emission line spectrum that can verify the tidal 
disruption scenario and constrain the mass of the BH. 
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1 INTRODUCTION 

With the exquisite angular resolution of Chandra, it is pos- 
sible to resolve low mass X-ray binaries (LMXBs) in nearby 
galaxies and to associate these LMXBs with globular clus- 
ters (GCs; e.g., Sarazin et al. 2000). This capability has en- 
abled the discovery of several black hole (BH) candidates 
in GCs outside of the Milky Way. Maccarone et al. (2007) 
presented the first strong candidate for a BH in a GC, an 
ultraluminous X-ray (ULX) source in the globular cluster 
RZ 2109 associated with the elliptical galaxy NGC 4472. 
The object's X-ray luminosity of Lx = 4 x 10 39 erg s _1 
is too large to be explained by accretion onto a neutron 
star (NS), and the source's strong variability rules out 
the possibility of a superposition of multiple NSs. Sub- 
sequently, three additional ULXs have been identified as 
GC BH candidates on the basis of their strong X-ray vari- 
ability; one in NGC 1399 (Shih et al. 2010), one in NGC 
3379 (Brassington et al. 2010), and a second candidate in 
NGC 4472 (Maccarone et al. 2011). Variable X-ray sources 
with luminosities above the Eddington limit of a NS of- 
fer the least ambiguous evidence of a BH accretor; how- 
ever, Barnard et al. (2011) argued for BH accretors in three 
LMXBs located in M31 GCs that exhibit low/hard-state 
spectra at luminosities near a NS's Eddington limit, and in 
a high luminosity recurrent transient in a fourth M31 GC. 

Irwin et al. (2010, hereafter 110) discussed an additional 
BH candidate in an NGC 1399 GC; however, the nature of 
this source is not as clear-cut as that of the other candidates. 
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The object's X-ray luminosity, Lx ~ 2 x 10 39 erg s _1 , sug- 
gests a BH accretor, but the 35% decline in the luminosity 
between 2000 and 2008 does not convincingly rule out a su- 
perposition of sources. The X-ray spectrum can be fitted 
with a power law of photon index F — 2.5, which is much 
softer than the slope of the power law used to fit the sum 
of low L x LMXBs (r ~ 1.6). This is also slightly softer 
than the spectra of other sources with Lx ~ 10 39 erg s _1 , 
which might indicate the presence of a more massive BH. 
Adding to the intrigue, optical spectroscopy of the NGC 
1399 GC revealed bright [N n] AA6548, 6583 and [O in] 
AA4959, 5007 emission lines with luminosities ~10 36 erg s _1 , 
but no Balmer emission lines. This is not the only GC BH 
with optical emission lines, Zepf et al. (2008) found bright, 
broad [O in] A5007 emission lines without accompanying 
Balmer lines in the optical spectrum of RZ 2109. 110 argued 
that these two sources could be evidence of the tidal dis- 
ruption or detonation of a white dwarf by an intermediate 
mass BH, based on the lack of hydrogen emission lines in 
the spectra. 

A BH in the core of a galaxy or a GC can tidally disrupt 
stars that pass too closely, and then accrete a portion of the 
debris (Hills 1975; Frank & Rees 1976; Sigurdsson & Rees 
1997; Baumgardt et al. 2004; Ramirez-Ruiz & Rosswog 
2009). For stars initially on parabolic orbits with the BH, 
roughly half of the mass of the disrupted star will become 
bound to the BH, and the other half will remain unbound 
and stream back out into the cluster (Lacy et al. 1982; 
Rees 1988). Roos (1992) first considered the emission lines 
produced when the unbound debris is photoionised by the 
UV/X-ray emission generated when the bound portion ac- 
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cretes onto the BH. Detailed studies of the emission lines 
produced when a main-sequence star is tidally disrupted 
by a supermassive BH have been made using numerical 
(Bogdanovic et al. 2004) and analytic (Strubbe & Quataert 
2009, 2011) models for the dynamical evolution of the debris. 
Clausen & Eracleous (2011) considered the case of a white 
dwarf being tidally disrupted by an intermediate mass BH 
and found that the emission lines observed in the two BH 
candidates described above were not consistent with those 
expected from a tidally disrupted white dwarf. 110 reported 
that the [N n] and [O in] doublets observed in the NGC 
1399 GC have roughly the same luminosity, but the mod- 
els predict that the [N n] doublet should be two orders of 
magnitude fainter than the [O in] doublet. The white dwarf 
tidal disruption models predict that the [O in] A5007 lumi- 
nosity will reach a peak value consistent with that observed 
in RZ 2109 a few years after tidal disruption. However, the 
X-ray source in RZ 2109 was detected with ROSAT in 1994 
(Colbert & Ptak 2002), which indicates that the observa- 
tions reported in Zepf et al. (2008) were made more than a 
decade after the BH began accreting. 

Other models have been proposed to explain the 
emission lines observed in GCs hosting BH candidates. 
Maccarone et al. (2010) found that the X-ray emission from 
RZ 2109 had varied by a factor of 20 over a period of six 
years, a timescale far too short to be consistent with a tidal 
disruption, and suggested that the source is a hierarchical 
triple system with an inner BH-white dwarf binary. Such 
a system is in line with the work of Ivanova et al. (2010), 
who suggested that a triple system may be necessary to 
form a BH-white dwarf binary. Ripamonti & Mapelli (2012) 
showed that nova ejecta that have been photoionised by a 
ULX could produce bright, broad [O m] emission lines sim- 
ilar to those observed in RZ 2109. Porter (2010) considered 
an accretion disc origin for the emission lines in RZ 2109 and 
the GC in NGC 1399. Maccarone & Warner (2011) argued 
that the lines observed in the NGC 1399 GC were formed in 
the photoionised winds of R Corona Borealis (RCB) stars. 
Here we propose that the X-ray continuum and optical emis- 
sion lines observed in the NGC 1399 are the result of a tidal 
disruption of a horizontal branch (HB) star. In section 2 we 
describe the ingredients of our models, in section 3 we de- 
scribe the results of our calculations, and in section 4 we 
compare our model with competing ones and discuss an ob- 
servational test to distinguish between them. 



2 MODELS 

We modelled the emission-line spectrum produced when 
a HB star on an unbound orbit is disrupted by a 
BH, using the method outlined in Clausen & Eracleous 
(2011), which is an adaptation of the procedure used by 
Strubbe & Quataert (2009). Our approach combines ana- 
lytic prescriptions for the dynamical evolution of the debris 
with Cloudy (Ferland et al. 1998) photoionisation models to 
compute the luminosities and line-profiles of the emission 
lines produced in the unbound portion of the debris. These 
models require as inputs the mass of the disrupted star M*, 
the radius of the disrupted star ii*, the mass of the BH Mbh, 
the pericentre distance of the orbit Tip , and the composition 
of the disrupted star. We determined the properties of the 
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Figure 1. Mass fractions (top) of carbon (solid), nitrogen (dot- 
ted), and oxygen (dashed) vs. time and radius (bottom) vs. time 
from the MESA star model of a 1 Mq star with Z = Zq/2. 

HB stars used in our models using the stellar evolution code 
MESA star (Paxton et al. 2011). 

2.1 Horizontal Branch Star Models 

A GCs horizontal branch morphology is strongly influenced 
by its metallicity, with the HBs in most metal rich GCs ly- 
ing redward of the instability strip and those of most low 
metallicity clusters lying to the blue side. However, it has 
been known for some time that metallicity alone cannot 
account for the range of HB morphologies observed, and 
that HB morphology is determined by metallicity and at 
least a second parameter (e.g., Sandage & Wallerstein 1960; 
van den Bergh 1967). By the mid 1990s, a number of stud- 
ies (e.g. Searle & Zinn 1978; Lee et al. 1994) had established 
cluster age as a likely solution to the second parameter prob- 
lem. Recent, independent studies by Dotter et al. (2010) and 
Gratton et al. (2010) have confirmed that the second param- 
eter is the age of the GC. The GC hosting the X-ray source 
considered here is quite red (B — I — 2.25; Kundu et al. 
2007), and therefore likely has a HB consisting primarily 
of red clump stars, which is typical of high metallically 
GCs. However, there is evidence that some metal rich GCs 
harbour extreme horizontal branch (EHB) and even blue 
hook stars (e.g. Dalessandro et al. 2008; Rey et al. 2007; 
Peacock et al. 2010), so we will consider the tidal disrup- 
tion of both red clump and EHB stars. 

Using MESA star, we evolved a 1 Mq star with Z = 
Zq/2 to determine the mass, radius, and composition of 
the red clump star. The star evolved to the horizontal 
branch after ~ 9.9 Gyr, at which point M* = 0.67 Mq 
and = 6.2 Rq. Furthermore, CNO cycle burning had 
enriched the core of the star with nitrogen at the expense 
of carbon and oxygen. Leading up to the HB phase, the to- 
tal mass fraction of nitrogen rose and those of oxygen and 
carbon fell as the star's relatively carbon and oxygen rich 
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Figure 2. Emission-line luminosity ratios over time from several 
tidal disruption models with varying Mgg. The long-dashed line 
in the top panel shows the constraint on the [N n]/[0 III] ratio 
measured by 110. Models that fall near this line are consistent with 
observations. In the lower two panels, the shaded region indicates 
line ratios ruled out by the observations of 110, only models that 
reach into the unshaded region are consistent with observations. 
For each model shown here, we considered a red clump star on 
an orbit with Rp = Rp. 



envelope was shed. The ratio of nitrogen to oxygen is high- 
est right after helium ignition, at the start of the HB phase. 
We used the mass, radius, and composition of the HB star 
at this point as our fiducial model. At this point the star 
was made up of 15% H, 84% He, 0.1% C, 0.5% N, and 
0.2% O, by mass, with all other elements present in their 
initial mass fraction. The mass and radius of our adopted 
HB star is consistent with previous models (e.g., Dorman 
1992; Charbonnel et al. 1996). Furthermore, measurements 
of the surface abundances in HB stars have revealed a ni- 
trogen enhancement, but the surface abundance of nitrogen, 
which depends on the stars mixing history, is highly variable 
from star to star (e.g., Behr et al. 1999; Gratton et al. 2000; 
Tautvaisiene et al. 2001). 

We found Ah = 0.49 M Q and = 0.12 Rq for an 



Figure 3. Emission-line luminosity ratios over time from several 
tidal disruption models with varying Mgn, tidal disruption pa- 
rameter /3, and stellar type. The long-dashed line in the top panel 
shows the constraint on the [N ll]/[0 III] ratio measured by 110. 
Models that fall near this line are consistent with observations. 
In the lower two panels, the shaded region indicates line ratios 
ruled out by the observations of 110, only models that reach into 
the unshaded region are consistent with observations. 

EHB star by running a MESA star model of a 1 Mq star with 
Z = 0.1 Zq. This combination of M* and R+ are in good 
agreement with the asteroseismology measurements of the 
masses and radii of hot subdwarfs (Charpinet et al. 2008; 
van Grootel et al. 2010). The composition of the EHB star 
was 5.2% H, 94.5% He, 5 x 10~ 3 % C, 0.1% N, and 0.03% O. 

2.2 Tidal Disruption Models 

Here we briefly describe our tidal disruption models, and 
we refer the reader to Clausen & Eracleous (2011) and 
Strubbe & Quataert (2009) for a complete description. We 
assume that the HB star is on a parabolic orbit with 
i?p < Rt, where Rt is the tidal disruption radius Rt = 
R*(Mbu/ R*) 1 ^ 3 . The strength of the encounter is set with 
the tidal disruption parameter /3 = Rt/Rp- After tidal dis- 
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ruption, half of the disrupted star falls back to the BH and 
forms a thin accretion disc. The rate at which material falls 
back declines with time as M oc t~ 5 ^ 3 . When M exceeds the 
BH's Eddington rate, we assume that the excess is driven 
away in an outflow but we do not include the outflowing 
material in our models (see Strubbe & Quataert 20f f for an 
exploration of the outflow.) This results in the accretion flare 
having a constant luminosity of L — 1.3 x 10 38 Mbh erg s _1 
during this super-Eddington phase. 

The other half of the debris remains unbound and 
moves away from the BH with a range of velocities given 
by v(4>) = (2GM B H/-Rp) 1/2 cot(0/2), where is the az- 
imuthal angle in the orbital plane and is in the range 
(tt - (12i?*/i? P ) 1/2 ) < <f> < 7T. The density in the debris 
tail is given by n((j>,t) = no[v(<j))t]~ 3 , where t is the time 
since tidal disruption and the normalisation factor no is set 
such that the total mass in the debris tail is M*/2. We di- 
vide the debris tail into six azimuthal zones and compute 
photoionisation models for each zone with version 10.00 of 
Cloudy (Ferland et al. 1998), using the density determined 
by the dynamical model and the composition given by the 
MESA star models. The spectral energy distribution (SED) 
and intensity of the ionising radiation are determined from 
the dynamical model of the bound debris. We assume a 
standard a-disc model with a = 0.1, an inner radius of 3 
Schwarzschild radii, an outer radius of 2 R p , and mass ac- 
cretion rate equal to the mass fallback rate described above. 
Using the range of temperatures in the accretion disc, we 
construct a multicolour blackbody which we us as the ion- 
ising continuum for the Cloudy models. We have not added 
an X-ray power law to the SEDs used in these models. Pho- 
toionisation by X-rays can result in significant heating of the 
unbound material and affect the luminosities of collisionally 
excited emission lines. However, Clausen & Eracleous (2011) 
found that the luminosities of the optical emission lines of 
interest here were not sensitive to this additional component. 

We modelled the emission lines produced with 20 differ- 
ent sets of initial conditions to explore the parameter space 
in BH mass (30 M < M B h < 1O 4 M ), HB star structure 
and composition (red clump and EHB), and tidal disruption 
parameter /3 (0.1 < /3 < 1). 



3 RESULTS 

Each combination of Mbh, HB star type, and p that we 
considered predicted [N n] A6583 and [O m] A5007 emission 
lines with luminosities ~ 10 36 erg s _1 at some point during 
the post tidal disruption evolution. However, only those with 
a red clump HB star, a BH with Mbh < 200 Mq , and f) ~ 1 
were able to reproduce the emission line luminosity ratios 
reported in 110. These constraints are [N n] A6583/Ha > 7, 
[O in] A5007/H/3 > 5, and [N n] A6583/[0 in] A5007~ 1. 
We plot the time evolution of these ratios for a subset of our 
models in Fig. 2. In the models plotted here, we have used 
our fiducial red clump HB star and assumed that Rp = Rt, 
but varied Mbh. 

In each of these models, the Balmer lines outshone the 
[O in] and [N n] lines initially. Then, as the debris cloud ex- 
panded and n decreased, the luminosity of the Balmer lines 
declined because the volume emissivity of these permitted 
transitions is proportional to the product of the electron and 
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Figure 4. FWHM of the [N II] A6583 emission line 125 years after 
tidal disruption, synthesised from a model using the red clump 
star, /3 = 1, and Mbh = 100 Mq. The angles do and io describe 
the orientation of the observer (see fig 1 in Clausen & Eracleous 
(2011)). The black curve marks the observed FWHM. Any pair of 
00 and io along the black line produces a FWHM that is consistent 
with observations. 



ionised hydrogen number densities. At the same time, the 
[O in] A5007 and [N n] A6583 lines became brighter. These 
forbidden lines are produced when ions are collisionally ex- 
cited into a metastable state and then de-excite radiatively. 
As the density in the unbound debris decreased, the colli- 
sion rate declined, allowing a larger fraction of the excited 
ions to decay radiatively, as opposed to collisionally, and the 
forbidden emission-lines grew brighter. As a result, all of the 
models shown in Fig. 2 can satisfy the [O m]/H/3 constraint, 
provided we are observing the debris more than 40 years 
after tidal disruption. The models with Mbh = 10 3 and 
10* Mq were unable to reproduce the observed [N n]/Ha 
ratio. These more massive BHs permit higher accretion rates 
that led to an increased flux of ionising photons. By the time 
the density in the debris was low enough to produce a signifi- 
cant [N ii] luminosity, [N n] emission was suppressed because 
the cloud had become so highly ionised that much of the N + 
was further ionised to N ++ . 

This over-ionisation of the unbound debris also drove 
the evolution of the [N n]/[0 in] ratio. In all of the models, 
the ratio steadily declined as an increasing portion of the 
oxygen and nitrogen in the cloud became doubly ionised. 
With Mbh = 10 3 or 10 4 Mq, the unbound debris became 
too highly ionised for the [N n]/[0 mjratio to be ~ 1 within 
ten years of tidal disruption (see Fig. 2). We found that in 
models with Mbh < 200Mq, the decline is slow enough 
that the luminosities of the two forbidden emission lines can 
differ by less than a factor of two for up to 200 years. 

Fig. 3, shows a number of alternative models for strong 
disruption and/or blue HB stars. These models fail to re- 
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Figure 5. Evolution of [N ll]/[0 III] for HB stars of different 
composition. We show two sets of models - young HB stars 30 
Myr after He ignition when the nitrogen and oxygen abundance 
are equal (M BH = 50Mq dashed; Af BH = IOOMq dotted) and 
the same models but with the HB star older so that 80 Myr 
of He burning has made oxygen twice as abundant as nitrogen 
(M BH = 50Mq dot-dashed; A/ BH = 100M Q solid). For a 100 
Mq BH, the observed [N n]/[0 III] ratio cannot be reproduced 
once helium burning has created an oxygen abundance that is 
twice the nitrogen abundance. 



produce the observed [N n]/[0 in] ratio for an extended pe- 
riod of time, in contrast to our models of weakly disrupted 
red HB stars. Like the models involving high mass BHs dis- 
cussed above, the unbound debris cloud becomes too highly 
ionised to produce significant [N n] emission. Using a simple 
estimate of the ionisation parameter U — %/njj, we can 
illustrate why some debris clouds become over ionised and 
others do not. Here, n 7 is the number density of photons 
with energies above 13.6 eV and uh is the hydrogen num- 
ber density. Recasting the ionisation parameter in terms of 

— 1/2 " *\ I A 

our model parameters, we find U oc i?* p M ' t. Debris 
from stars with a larger kinetic energy at pericentre pas- 
sage, due to a more compact star with smaller Rt, an orbit 
with ft > 1, or a more massive BH, spanned a wider range 
in velocity and quickly expanded to a low density, resulting 
in a large U. As mentioned above, more massive BHs lead 
to a higher flux of ionising photons, also contributing to in- 
creased U. During the super-Eddington phase, M oc M B h 
was constant and the ionisation parameter increased with 
time. The N + fraction reached a maximum 1 when U ~ 0.06, 
and at larger values of U, the luminosity of the [N n] line 
declined. Keeping the ionisation state in the unbound de- 
bris low enough to produce significant [N n] emission after 
the super-Eddington phase had ended required the relatively 
large radius of the red clump star and a BH of relatively low 
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Figure 6. Light curves from a model using the fiducial red clump 
HB star, /3 = 1, and M B jj = 100 Mq. Emission line light curves 
are shown for N V A1240 (grey, filled circle), C IV A1548 (green, 
plus), [O III] A5007 (purple, asterisk), [N II] A6583 (light blue, 
triangle), N II A5755 (red, filled box), O III A4363 (dark blue, dia- 
mond), Ha (black, open box), H/3 (orange, x). The dashed curve 
shows Lx scaled down by a factor of 2000. Note the bright N V 
and C IV doublets predicted by the model. The dotted portion 
of each curve shows the line luminosity while the fallback rate is 
super-Eddington and our model is uncertain. 

mass, Af B H < 200 Mq. After the super-Eddington phase, U 
declined as t~ 1//4 and N ++ recombined into N + . This drove 
the [N n]/[0 in] ratio back towards unity in some cases (see 
models with /3 > 1 in Fig. 3). However, at this point the 
density in the cloud is so low that the luminosities of the 
lines are well below the observed values. 

Finally, by combining the dynamical model with the 
photoionisation model, we synthesised the emission line pro- 
files of the [N n] and [O in] lines. Fig. 4 shows the range 
of observed full width at half maximum (FWHM) of the 
[N n] line produced in a the model of a red clump star on 
an orbit with (3 = 1 being tidally disrupted by a IOOA/q 
BH, 125 years after tidal disruption. The angles 6q and io 
describe the orientation of the observer, and are the angle 
between the line of site and the direction of pericentre and 
the orbital plane, respectively. The black line shows the ob- 
served FWHM = 140 km s _1 . For a range of observer ori- 
entations, the FWHM of the synthesised emission lines is 
consistent with the observations. 



The value of U at which this maximum occurs is abundance 
dependent. The value given is for the red clump star composi- 
tion. The maximum occurs at a much smaller value of U if solar 
abundances are used. 
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4 DISCUSSION 

We have shown that the photoionised debris of a tidally dis- 
rupted HB star can account for the emission lines observed in 
optical spectra presented in 110. Reproducing the line ratios 
within the context of this model requires Mbh ^200 Mq. 
In these scenarios, the super-Eddington phase lasts between 
90 and 155 years. Since we have not included the outflow 
produced during this phase in our models, the emission 
line luminosities calculated during this period are uncertain. 
However, the modelled emission line luminosities remain 
consistent with the observations after the super-Eddington 
phase has ended. Furthermore, the X-ray luminosity of the 
accretion flare does not decline to the observed value of 
Lx ~ 2x 10 39 erg s _1 until after the mass fall back rate has 
dropped below the Eddington rate. Archival ROSAT HRI 
images establish that the X-ray source has been on since at 
least Jan. 1996. At times greater than t — 200 years after 
disruption, our assumption of steady state accretion begins 
to break down because the viscous time in the accretion disc 
approaches t. We have also checked that the conditions in 
the photoionised debris change slowly enough that the un- 
bound debris remained in photoionisation equilibrium for 
the entire period explored in our models. 

We did not consider spinning BHs in this study. The 
last stable orbit around a spinning BH is closer to the event 
horizon, allowing the accretion disc to extend to smaller radii 
than the discs considered here. Accretion discs around spin- 
ning BHs produce a larger flux of ionising photons with an 
SED that peaks at a higher photon energy because the max- 
imum temperature in these discs is larger than that of ac- 
cretion discs around non-spinning BHs. However, unlike the 
supermassive BHs found in the cores of galaxies, BHs in GCs 
have not gained an appreciable portion of their mass though 
disc accretion. Belczynski et al. (2008) found that any accre- 
tion that does occur while the BH is in a mass transferring 
binary does not substantially increase its spin. However, it 
is possible that the GC BHs are born with a large spin. If 
this is the case, the HB tidal disruption scenario remains 
a plausible explanation for the observed emission lines, but 
the BH mass inferred from the models would change by a 
factor of a few. 

The [N n]/[0 m] ratio is sensitive to the abundance of 
nitrogen and oxygen in the HB star. In the models discussed 
above, we used the composition of the HB star at the onset 
of helium burning, when the ratio of nitrogen to oxygen is 
largest. Helium burning destroys nitrogen and creates oxy- 
gen, so N/O declines as the star evolves on the HB. Fig. 5 
shows how [N n]/[0 in] responds to changes in composition 
of the HB star. In these models, we used the composition 
of the red clump star at later times, as computed in the 
MESA star models, corresponding to points where the nitro- 
gen and oxygen abundance are equal, and where the nitro- 
gen abundance is half that of oxygen. The mass and radius 
of the star did not change appreciably in this time, so we 
used the same values as in our fiducial model; and we used 
13 — 1. When Mbh = 50 Mq, the models with both abun- 
dance sets are consistent with the observed [N n]/[0 in]. The 
[N n]/[0 in] ratio drops too quickly to be consistent with ob- 
servations when Mbh = 100 Mq and the oxygen abundance 
is twice that of nitrogen. This sets the period during which 
the HB star can be disrupted by the BH and produce the 



observed emission line spectrum. For the 100 Mq BH, the 
HB star must be disrupted within 30 Myr of helium ignition 
to satisfy the [N n]/[0 in] constraint, while a 50 Mq BH 
can disrupt the HB star within 80 Myr of helium ignition 
and meet the constraint. 

It is appealing that the tidal disruption models require 
a red horizontal branch star to meet the observational con- 
straints on the emission line ratios. The host globular cluster 
is red and likely exhibits the red horizontal branch mor- 
phology typical of a metal rich cluster. However, given the 
predicted rate at which GC BHs tidally disrupt stars, it is 
surprising that such an event might have been observed. 
Employing the tidal disruption rate equation computed by 
Baumgardt et al. (2004), with the mass and radius of a 
red clump HB star, Mbh = 100 Mq, a core density of 
5 x 10 5 pc~ 3 , a core velocity of 10 km s -1 , and assuming 
that 5% of the stars in a globular cluster's core are HB stars, 
we find a HB tidal disruption rate of 1.2 x 10" 10 yr _1 GCT 1 . 
Given that the bright [N n] and [O in] emission lines persist 
for 200 years, and that the space density of globular clus- 
ters is ~ 6 Mpc~ 3 (Brodie & Strader 2006; Croton et al. 
2005), at any time there should only be 5 x 10~ 3 observable 
HB tidal disruptions within the 20 Mpc distance to NGC 
1399. Furthermore, if the NGC 1399 GC source is the af- 
termath of a HB tidal disruption, it is also puzzling that 
several main sequence tidal disruptions have not been ob- 
served because they because they should occur more fre- 
quently (Baumgardt et al. 2004). 

The low disruption rate and the dearth of main sequence 
tidal disruptions are only an issue if the tidal disruption oc- 
curs through the conventional loss cone orbit channel. One 
possibility that would result in the preferential disruption 
of post main sequence stars is for the disrupted star to be 
bound to the BH. If this is the case, tidal disruption could be 
triggered in two different scenarios; 1) a scattering induced 
merger or 2) a triple decay. In a scattering induced merger, 
the HB star progenitor is bound to the BH, expands to fill 
its Roche Lobe while on the giant branch, and transfers ma- 
terial to the BH. During mass transfer, the semi-major axis 
of the binary expands, greatly increasing the likelihood that 
the binary will undergo an encounter with another star in 
the GC. Such an encounter could impulsively increase the 
eccentricity and decrease the semi-major axis of the binary 
and result in the tidal disruption of the bound star (Davies 
1995; Sigurdsson & Phinney 1993, 1995). Successive distant 
encounters can also to drive the eccentricity high enough for 
Rp = Rt- We have estimated the timescale for this process 
using the expression given in Maccarone (2005), which uses 
the encounter cross sections computed by Heggie & Rasio 
(1996). In a binary with a post-mass-transfer semi-major 
axis of 1 AU containing a 100 Mq BH, a 0.65 Mq HB 
star will be driven to the tidal disruption radius in 10 Myr. 
This is well within the 30 Myr window required to produce 
[N n] and [O in] emission lines of similar luminosity. Here 
we have assumed that the field stars have mass 0.8 Mq, a 
cluster density of 5 x 10 s pc -3 , and a velocity of dispersion 
of 10 km s _1 . A scattering induced merger would be more 
likely in a dense GC. 

In the triple decay, the BH and HB star progenitor are 
initially the inner binary in a stable, hierarchal triple sys- 
tem. Again, evolution of the HB progenitor results in a phase 
of mass transfer and expansion of the inner binary's semi- 
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major axis. In this case, expansion of the inner binary drives 
the triple to instability and triggers the tidal disruption. A 
similar scenario was discussed by Perets & Kratter (2012), 
only these authors invoked mass loss, not mass transfer, to 
drive the expansion of the inner binary. The triple decay 
scenario would require a medium density GC, because triple 
systems are unstable to external perturbations in dense GCs. 
Both the scattering induced merger and triple decay scenar- 
ios increase the odds that tidal disruptions are linked to post 
main sequence stars, but determining if disruption through 
either scenario would occur more frequently than the loss 
cone tidal disruptions discussed above is beyond the scope 
of this paper. 

Our models also predict several strong UV emission 
lines. These lines could be used to test the tidally disrupted 
HB star hypothesis for the ULX in NGC 1399 and to fur- 
ther constrain the mass of the BH. The strongest lines are 
N v AA1239, 1243, O vi AA1032, 1038, N iv] A1486, and C iv 
AA1548, 1550. These lines should have FWHM similar to the 
observed optical lines. The bright C iv doublet could be used 
to constrain the composition of the disrupted star and nar- 
row the allowed range in Mbh. The light curves for the N v 
and C iv emission lines are shown in Fig. 6 along with the 
emission line light curves for the optical lines with observa- 
tional constraints. In contrast, Maccarone & Warner (2011) 
do not report any strong UV lines from their models of pho- 
toionised RCB star winds. 



5 CONCLUSION 

We have presented models for the emission lines generated in 
the photoionised debris of tidally disrupted HB stars. Within 
the assumptions made in our models, the emission line spec- 
tra produced when a red clump star is disrupted by a BH 
with Mbh < 200 are consistent with the optical emission 
lines observed by 110 in the ULX hosting globular cluster 
in NGC 1399. The models also make testable predictions 
about the source's UV-emission line spectrum. Observations 
of the UV-emission lines could significantly improve our un- 
derstanding of this source and are likely possible with the 
Hubble Space Telescope. 

ACKNOWLEDGMENTS 

We thank the anonymous referee for a prompt and helpful 
review. DC is supported by the PSU Academic Computing 
Fellowship. The authors thank the Aspen Center for Physics 
for its hospitality. 



REFERENCES 

Barnard R., Garcia M., Li Z., Primini F., Murray S. S., 

2011, ApJ, 734, 79 
Baumgardt H., Makino J., Ebisuzaki T., 2004, ApJ, 613, 

1143 

Behr B. B., Cohen J. G., McCarthy J. K., Djorgovski S. G., 

1999, ApJL, 517, L135 
Belczynski K., Taam R. E., Rantsiou E., van der Sluys M., 

2008, ApJ, 682, 474 
Bogdanovic T., Eracleous M., Mahadevan S., Sigurdsson 

S., Laguna P., 2004, ApJ, 610, 707 



Brassington N. J. et al., 2010, ApJ, 725, 1805 
Brodie J. P., Strader J., 2006, ARA&A, 44, 193 
Charbonnel C, Meynet G., Maeder A., Schaerer D., 1996, 

A&AS, 115, 339 
Charpinet S., van Grootel V., Reese D., Fontaine G., Green 

E. M., Brassard P., Chayer P., 2008, A&A, 489, 377 
Clausen D., Eracleous M., 2011, ApJ, 726, 34 
Colbert E. J. M., Ptak A. F., 2002, ApJS, 143, 25 
Croton D. J. et al., 2005, MNRAS, 356, 1155 
Dalessandro E., Lanzoni B., Ferraro F. R., Rood R. T., 

Milone A., Piotto G., Valenti E., 2008, ApJ, 677, 1069 
Davies M. B., 1995, MNRAS, 276, 887 
Dorman B., 1992, ApJS, 81, 221 
Dotter A. et al., 2010, ApJ, 708, 698 

Ferland G. J., Korista K. T., Verner D. A., Ferguson J. W., 
Kingdon J. B., Verner E. M., 1998, PASP, 110, 761 

Frank J., Rees M. J., 1976, MNRAS, 176, 633 

Gratton R. G., Carretta E., Bragaglia A., Lucatello S., 
D'Orazi V., 2010, A&A, 517, A81 

Gratton R. G., Sneden C, Carretta E., Bragaglia A., 2000, 
A&A, 354, 169 

Heggie D. C, Rasio F. A., 1996, MNRAS, 282, 1064 

Hills J. G., 1975, Nat, 254, 295 

Irwin J. A., Brink T. G., Bregman J. N., Roberts T. P., 

2010, ApJL, 712, LI 
Ivanova N., Chaichenets S., Fregeau J., Heinke C. O., Lom- 

bardi, Jr. J. C, Woods T. E., 2010, ApJ, 717, 948 
Kundu A., Maccarone T. J., Zepf S. E., 2007, ApJ, 662, 

525 

Lacy J. H., Townes C. H., Hollenbach D. J., 1982, ApJ, 
262, 120 

Lee Y.-W., Demarque P., Zinn R., 1994, ApJ, 423, 248 
Maccarone T. J., 2005, MNRAS, 364, 971 
Maccarone T. J., Kundu A., Zepf S. E., Rhode K. L., 2007, 
Nat, 445, 183 

Maccarone T. J., Kundu A., Zepf S. E., Rhode K. L., 2010, 

MNRAS, 409, L84 
Maccarone T. J., Kundu A., Zepf S. E., Rhode K. L., 2011, 

MNRAS, 410, 1655 
Maccarone T. J., Warner B., 2011, MNRAS, 410, L32 
Paxton B., Bildsten L., Dotter A., Herwig F., Lesaffre P., 

Timmes F., 2011, ApJS, 192, 3 
Peacock M. B., Maccarone T. J., Knigge C, Kundu A., 

Waters C. Z., Zepf S. E., Zurek D. R., 2010, MNRAS, 

402, 803 

Perets H. B., Kratter K. M., 2012, arXiv (1203.2914) 
Porter R. L., 2010, MNRAS, 407, L59 
Ramirez-Ruiz E., Rosswog S., 2009, ApJL, 697, L77 
Rees M. J., 1988, Nat, 333, 523 
Rey S. C. et al., 2007, ApJS, 173, 643 

Ripamonti E., Mapelli M., 2012, MNRAS, in press 

(arXiv: 1203.4233) 
Roos N., 1992, ApJ, 385, 108 
Sandage A., Wallerstein G., 1960, ApJ, 131, 598 
Sarazin C. L., Irwin J. A., Bregman J. N., 2000, ApJL, 544, 

L101 

Searle L., Zinn R., 1978, ApJ, 225, 357 

Shih I. C, Kundu A., Maccarone T. J., Zepf S. E., Joseph 

T. D., 2010, ApJ, 721, 323 
Sigurdsson S., Phinney E. S., 1993, ApJ, 415, 631 
Sigurdsson S., Phinney E. S., 1995, ApJ, 99, 609 
Sigurdsson S., Rees M. J., 1997, MNRAS, 284, 318 



© 0000 RAS, MNRAS 000, 000-000 



8 Clausen et al. 



Strubbe L. E., Quataert E., 2009, MNRAS, 400, 2070 
Strubbe L. E., Quataert E., 2011, MNRAS, 415, 168 
Tautvaisiene G., Edvardsson B., Tuominen I., Ilyin I., 2001, 

A&A, 380, 578 
van den Bergh S., 1967, AJ, 72, 70 

van Grootel V., Charpinet S., Fontaine G., Green E. M., 
Brassard P., 2010, A&A, 524, A63 
Zepf S. E. et al., 2008, ApJL, 683, L139 



© 0000 RAS, MNRAS 000, 000-000 



